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Abstract 

In this paper we study a novel spin chain with nearest-neighbors interactions de- 
pending on the sites coordinates, which in some sense is intermediate between the 
Heisenberg chain and the spin chains of Haldane-Shastry type. We show that when 
the number of spins is sufficiently large both the density of sites and the strength of 
the interaction between consecutive spins follow the Gaussian law. We develop an 
extension of the standard freezing trick argument that enables us to exactly com- 
pute a certain number of eigenvalues and their corresponding eigenf unctions. The 
eigenvalues thus computed are all integers, and in fact our numerical studies evi- 
dence that these are the only integer eigenvalues of the chain under consideration. 
This fact suggests that this chain can be regarded as a finite-dimensional analog of 
the class of quasi-exactly solvable Schrodinger operators, which has been extensively 
studied in the last two decades. We have applied the method of moments to study 
some statistical properties of the chain's spectrum, showing in particular that the 
density of eigenvalues follows a Wigner-like law. Finally, we emphasize that, unlike 
the original freezing trick, the extension thereof developed in this paper can be ap- 
plied to spin chains whose associated dynamical spin model is only quasi-exactly 
solvable. 
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1 Introduction 



Solvable spin chains have enjoyed a growing popularity in the last few years, 
due in part to their novel applications to SUSY Yang-Mills and string theo- 
ries [Tp|3Pll5]|6] . The prime example of such chains is the celebrated Heisenberg 
model describing spins in a one-dimensional lattice with isotropic nearest- 
neighbors interactions independent of the site. The Hamiltonian of the model 
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is given by [7] 

i 

where = {Sf^S'^^S^) is the spin operator of the i-th site, the sum runs 
from 1 to (as always hereafter), and Sn+i = Si. As is well-known, for 
spin 1/2 the model ([T]) can be exactly solved using Bethe's ansatz pi?|ITU] . 
Several (partially) solvable generalizations of the Heisenberg chain ([T]) with 
short range interactions (at most between next to nearest neighbors) have 
been subsequently proposed in the literature. These include, in particular, 
the family of chains with arbitrary spin and nearest-neighbors interactions 
polynomial in Sj ■ Sj+i of Refs. pTfT^ , as well as several models whose ground 
state can be written in terms of "valence bonds" [TSffH] . 

A different type of solvable spin chain with long-range position-dependent 
couplings was introduced independently by Haldane and Shastry 
This chain describes a system of A^ spins equally spaced on a circle, such 
that the strength of the interaction between each pair of spins is inversely 
proportional to their chord distance. The motivation for introducing the HS 
chain was the fact that its exact ground state coincides with Gutzwiller's 
variational wave function for the Hubbard model [T7lll8fl9j when the strength 
of the on-site interaction tends to infinity. We shall write the Haldane-Shastry 
(HS) Hamiltonian as 

1 ?7r 

HhS = 2 E - - S^J) , = IT^ ' (2) 

i<j 

where Sij is the operator exchanging the i-th and j-th spins. Although the 
particles' spin in the original HS chain was assumed to be 1/2, one can more 
generally consider particles with n internal degrees of freedom transforming 
under the fundamental representation of su(n). In this case the spin permu- 
tation operators can be written in terms of the fundamental su(n) generators 
at each site k (normalized such that tr( jf ) = ^ S"'^) as 

a=l 

Note that for spin 1/2 particles {n = 2), we have Sj = (J/, Jf, Jf). 

The HS chain is naturally related to the scalar Sutherland model of An type 
[20|[2T] and its spin version introduced in Refs. [22ll23|24j . In fact, Polychron- 
akos [25] noted that the complete integrability of the HS chain could be de- 
duced from that of the spin Sutherland model by suitably taking the strong 
coupling limit (the so-called "freezing trick"). Moreover, the latter author 
applied this technique to construct an integrable spin chain related to the 
Calogero (rational) model of type [26j. The Hamiltonian of this chain, 
usually referred to in the literature as the Polychronakos-Frahm (PF) chain. 
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reads 



^pF = E(0-0)-'(i-5.,), (3) 

i<j 

where the chain sites Q are now the equihbrium positions of particles in 
the scalar part of the potential of the Calogero spin model of type. In 
particular, the sites of the PF chain are not equally spaced, unlike those of 
the HS chain. Indeed, Frahm [27] pointed out that the sites of the PF chain 
are the zeros of the iV-th Hermite polynomial, which satisfy the system of 
algebraic equations 

= E7^' z = l,...,iV. (4) 



It turns out that both the HS and the PF chains, featuring long-range position- 
dependent interactions, can be solved in a more detailed and explicit way 
than the chains of Heisenberg type, characterized by the short range and 
position independence of the interactions. For instance, the partition function 
of the models ([2]) and ([3]) can be evaluated in closed form for arbitrary 
using the freezing trick [28,29], and its expression is relatively simple in both 
cases. The spectrum, which consists of a set of integers (consecutive in the 
case of the PF chain), is highly degenerate due to an underlying Yangian 
symmetry [30] . Moreover, it has been recently shown that even for moderately 
large N the density of eigenvalues of the HS chain is Gaussian to a high 
degree of approximation, and that the density of spacings between consecutive 
levels follows a simple distribution different from the usual Poisson or Wigner 
laws [221 ■ In fact, there is strong evidence that these results also hold for other 
spin chains of Haldane-Shastry type, as e.g. the BCn or the supersymmetric 
versions of the original HS chain [3T|32] . 

In this paper we shall consider a novel type of spin chain, which in some sense 
is intermediate between the Heisenberg and the Polychronakos-Frahm chains. 
Its Hamiltonian is obtained from that of the PF chain Q-dl]) by retaining 
only nearest-neighbors interactions, namely 

^ = Ete-em)~Ml-^M+i), (5) 

i 

where the sites are defined by restricting the sum in Eq. (jlj) to nearest 
neighbors: 

+ , , t = l,...,N. (6) 

In the previous equations we are identifying S^^^^i with Sm, ^n+i with 
and ^0 with ^^v- 
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We shall see that the chain possesses several remarkable properties, whose 
study is the purpose of this paper. In the first place, the spin chain ([5]) is 
related along the lines of the freezing trick to the spin dynamical model 



+ T.7Zr^Z-^i^-S^,^+l), (7) 

where = J^i^h ^ > 1/2, and we have identified xq = xn and xn+i = xi. 
We have shown in our recent papers [33f3^ that an infinite proper subset of 
the spectrum of the Hamiltonian ([7]) can be computed in closed form, so that 
this model is quasi-exactly solvable (QES) [35|36II37] . The Hamiltonian ([7]) is 
a spin version of the QES scalar model 

Hsc = i/|5,,,+i^l (8) 

introduced in [38] by Auberson, Jain and Khare. In particular, it was shown 
in the latter reference that 

;i = e-t^'[||xi -Xi+il" (9) 

i 

is the ground state function of the model ([8]), with eigenvalue Eq = Na{2a+1). 
We shall prove that the sites of the chain are in fact the coordinates of 
the unique maximum of /i in the domain 

C=[xeR^ \ xi<---<xn}. (10) 

By numerically solving Eq. we shall see that for sufficiently large the 
chain sites are normally distributed with zero mean and unit variance. We 
shall present a simple deduction of this property based on the analysis of the 
continuous limit of the algebraic system ([6]) as — > cxo. As a byproduct, we 
shall obtain an analytic formula providing a very accurate approximation to 
the sites' coordinates, valid even for moderately large values of A^. Another 
nontrivial consequence of this formula is the fact that for large A^ the strength 
of the coupling between the spins i and i + 1 as a function of their mean 
coordinate (^i + ^j+i)/2 also follows the Gaussian law, but with zero mean and 
variance 1/2. 

The spectral properties of the spin chain Hamiltonian are also remarkable. 
Indeed, by a suitable modification of the freezing trick one can show that Ti. 
possesses the eigenvalues 0,1,2 for arbitrary values of A^ and n, and exactly 
compute their corresponding eigenstates. Our numerical simulations evidence 
that these energies are the three lowest ones, and that for spin 1/2 none of 
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the remaining eigenvalues of Ti are integer^]. For n > 2, the spectrum of 
Ti also contains the integer eigenvalue 3, which appears to be not the fourth 
but the fifth lowest energy. The above properties suggest that the model ^ 
could be regarded as a quasi-exactly solvable chain, in the sense that only 
a certain number of eigenvalues and their corresponding eigenvectors of the 
Hamiltonian Ti can be computed in closed form. 

We have also studied the distribution of energy levels of the chain ^ for 
a large number of particles. Since the partition function of this chain is not 
known, we have performed a numerical calculation of the density of levels 
using the methods of moments [39P0] . It turns out that, in contrast with the 
typical behavior of spin chains of Haldane-Shastry type, the distribution of 
levels clearly deviates from the Gaussian law. As an indication of the accuracy 
of the approximate level density derived via the moments method, we have 
compared its mean and variance with the exact values obtained by taking 
traces of suitable powers of the Hamiltonian. From this discussion it also 
follows that for large the mean and variance of the energy behave as and 
, respectively, just as for the trigonometric chains of HS type in Refs. [2SEI]- 

The paper is organized as follows. In Section [2] we study the distributions 
of the chain sites and the couplings, comparing the results obtained with 
those for the PF chain. Section [31 which is the core of the paper, is devoted 
to the determination of the integer eigenvalues of the chain ^ and their 
corresponding eigenstates. We also present in this section a detailed example 
for the case of 5 particles of spin 1, which motivates a number of conjectures on 
the degeneracy of the integer levels. In SectionlHwe use the method of moments 
to approximately compute the density of levels of the chain 7i, showing that it 
follows a Wigner-like law. Finally, in Section O we summarize our conclusions 
and outline possible future developments. For the reader's convenience, in 
Appendix [A] we present some background material on the exact eigenfunctions 
of the spin dynamical model ([7]) used in Section [3], while in Appendix [B] we 
include an overview of the method of moments. 



2 The chain sites 

We shall start this section by proving that the sites of the chain ([5]) are the 
coordinates of the unique maximum in the domain ffTOj) of the ground state 



^ As a matter of fact, the previous assertion does not hold for the cases = 3 
(for which the chain ([S]) reduces to the PF chain, whose eigenvalues are known to 
be integers), and = 4 (for which our numerical simulations indicate that all the 
eigenvalues are also integers). Therefore, in the rest of the paper we shall exclude 
these special cases from our discussion. 
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function (jO]) of the scalar Hamiltonian (jS]). It is convenient to write the ground 

state as 



/i 



where 



A(x) = ^log 



^i+l I 



has the same extrema as fi and is independent of a. Thus the equations 
defining the chain sites are just the conditions for ^ = (^i, . . . to be a 
critical point of A. The existence of a maximum of A in C is clear, since it 
is continuous in C and tends to — cxd both on its boundary and as r — > oo. 
Uniqueness follows from the fact that the Hessian of A is negative definite 
in C. Indeed, by Gerschgorin's theorem [41, 15.814], the eigenvalues of the 
Hessian of A at x lie in the union of the intervals 



dxf 



li 



dx] 



+ 7i 



where 7i = ^ 



dxidxj 



Since 



S^A 



Xi+l) ^ — {Xi — ^, 

d'^X 



dxidxi±i 



[^Xi Xi-^\j 



dxidxj 



if j^t,t±l, 



we have 



dx'} ^* 



-1 



and thus all the eigenvalues of the Hessian of A are strictly negative. 
Since A and C are invariant under the transformation 

Xi^^-XN-i+i, i = l,...,N, 
and A has a unique maximum in C, it follows that 

6 = —^N-i+l, 



:iii 



so that the chain sites are symmetric about the origin. In particular, the center 
of mass of the spins vanishes, i.e.. 



(12) 



We have numerically solved equations ([6]) for the chain sites for up to 250 spins. 
Before presenting our conclusions, a remark on the configuration of these sites 
is in order. The attentive reader might have been surprised by our claim that 
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the chain ([5]) features only nearest-neighbors interactions, in spite of the fact 
that the first spin interacts with the last one. We can avoid this objection by 
regarding the site coordinate as an arc length in a circle of radius 2^]y/ii, see 
Fig.dJ In this way the spins at the sites and C,n are indeed nearest-neighbors 
and, moreover, the strength of the interactions are inversely proportional to 
the squared distance between consecutive spins, measured along the arc. 




Fig. 1. Sites of the chain TC for = 20 spins. 



It is apparent from Fig. [T] (and also follows immediately from Eq. ([6])) that the 
sites C,i are not equally spaced. In fact, our computations show that for large 
values of the sites C,i follow with great accuracy a Gaussian distribution 
with zero mean and unit variance. More exactly, the cumulative density of 
sites (normalized to unity) 

^(a;)=iV-i5]^(x-6), (13) 

i 

where 6 is Heaviside's step function, is approximately given by 

1 



1 + erf(x/v^) 



(14) 



The agreement between the functions !F and F is remarkably good for > 100 
(see Fig. [2]for the case N = 150) and increases steadily with A^, e.g., the mean 
square error of the fit for 100, 150 and 200 spins are respectively 2.6 x 10~^, 
1.1 X 10-5 and 7.9 x 10"^. 

The fact that for large A^ the cumulative density of sites is well approximated 
by the Gaussian law (fT4|) can be justified by the following heuristic argument. 
Let x{t,N) be a smooth function such that x{i,N) = for i = 1,...,N, 
and define the rescaled function y{s,e) = x{s/e, 1/e). By Eq. IQ, the latter 
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T{x),F{x 




Fig. 2. Cumulative distribution functions J-^{x) (at its discontinuity points) and 
F{x) (continuous grey line) for N = 150 spins. 



function must satisfy the relation 

1 



y{s,e) -y{s - e,e) y{s,e) ~ y{s + e,e) 



yis,e) 



(15) 



for e = 1/A^ ^ 1 and s = j^, . . . ,1. Let us now assume that Eq. (fT5|) holds 
for all s G M and all e <^ 1. Writing 



k=0 



and using the expansion 

y{s,e) -y{s±e,e) = TZ/o(s)e 



±y[is))e' + 0{e' 



the leading term in Eq. ( |T5i) yields the differential equation 



// 1 2 

yo = yo yo 



The general solution of this equation is implicitly given by 

s = cq + Ci erf {yo{s)/V2^ . 



(16) 



Hence, up to terms of order e = 1/A^, the cumulative distribution function 
of the chain sites (normalized to unity) is approximated by the continuous 
function 

F{x) = Co + ci erf (x/ V2) . 

The normalization conditions F{—oo) = and F{oo) = 1 imply that Cq = 
ci = 1/2, and thus the empiric law (IT^ is recovered. 
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From Eq. (fT6l) (with cq = Ci = 1/2) it follows that the site can be determined 
up to terms of order 1/N hj the formula 



~ V2 erf 



_i /2A; - 
V N 



(17) 



If the sites were exactly given by the previous formula, they would satisfy 
the identity 

which is clearly inconsistent with the exact relation (ITTl) . However, the slightly 
modified formula 



~ erf 



-1 



2k -N 



(18) 



differs from ( JTTl) by a term of order 1/A^ and is fully consistent with the 
relation ( ITT]) . Although both ( fT7|) and ( |T8l) provide an excellent approximation 
to the chain sites for large A^, the latter equation is always more accurate than 
the former, and can be used to estimate C,k with remarkable precision even for 
relatively low values of A^, cf. Fig. [31 



Fig. 3. Sites coordinates and their continuous approximation (jlSp for A^ = 20 
spins. 



It is also of interest to determine whether the position of the last spin tends 
to infinity as A^ ^ 00, since according to our interpretation of the chain's 
geometry the number 2C,n/it is the radius of the circle on which the spins 
lie. From Eq. f[T8l) it follows that for large A^ the last spin's coordinate C,n is 
approximately given by 

^N^V2 err' (1 - ^) , (19) 

so that should diverge as A^ — 00. Of course, this assertion should be taken 
with some caution, since in Eq. (1181) the argument of the inverse error function 
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is correct only up to terms of order In order to check the correctness of 

the approximation (fT9|) . we recall the asymptotic expansion of erf~^('u) for 
M ^ 1 in Ref. [42j to replace fll9l) by the simpler formula 



N 



2ri — log 7] . 



(20) 



where 



Tj = lo^ 



As can be seen in Fig. HI the approximate formula fl20|) qualitatively reproduces 
the growth of when N ranges from 100 to 250. A greater accuracy can 
be achieved by introducing an adjustable parameter in Eq. f[T^ through the 
replacement 1/N ^ a/N, so that in Eq. (12(1 rj becomes 



T] = log 



N 



a\/n 



(21) 



In Fig. m we have also plotted the law (I2UI) - (I^ with the optimal value a 
0.94, which is in excellent agreement with the numerical values of for 
100,105, ...,250. 



N 

^wO 120 140 IGO 180 200 220 240 

Fig. 4. Position of the last spin ^at for = 100, 105, . . . , 250 and its continuous 
approximation ([20|)-([2T|) for a = 0.94 (solid line) and a = 1 (dashed line). 



The last property of the spin chain that we shall analyze in this section 
is the dependence of the coupling between neighboring spins on their mean 
coordinate. Calling 



(a-e 



we shall now see that when > 100 the Gaussian law 

hu ~ — e ^fc 
27r 



(22) 



(23) 
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holds with remarkable precision, cf. Fig. [51 Indeed, if x = x{k) denotes the 
RHS of Eq. ([18]) we have 

2k = Neif (^^^ +N+1, 



so that 



dx v27r i,.2 
e2 



dA; 

is of order 1/A^. Hence, up to terms of order 1/A^ we have 



(24) 



^ 2/ ~^V^ 2 



-2 



dx 



dk 



-2 



iV2 



27r 



(25) 



Since (up to terms of order l/N"^ 



x{k) - - 



Eq. (!23|) follows from ([25!) replacing /c by A; + |. 
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Fig. 5. Plot of the coupling between the spins k and k + 1 versus their mean 
position, fitted by the Gaussian ([23]) . for = 100 spins. 



We shall finish this section with a brief comparison of the previous properties 
with those of the PF chain ([3]). For large A^, the density of sites of the PF 
chain (normalized to unity), that is the density of zeros of the iV-th Hermite 
polynomial, is asymptotically given by the circular law [13] 

Pn{x) = :^V2N-x^. (26) 

The last site Ctv of the PF chain grows with N much faster than the corre- 
sponding site C,N of the chain ([5]), for the largest zero of the iV-th Hermite 
polynomial behaves as \/2N + 0{N'^^^); see, e.g., Finally, in contrast 
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with the PF chain, the polynomials determined for each N ^ 2 by the corre- 
sponding sites of the chain ([5]) do not form an orthogonal family. In other 
words, the polynomials 

pW(z)=n(^-6) (27) 

i 

do not satisfy a three-term recursion relation of the form 
as it can be verified using the explicit expressions 



3 Quasi-exact solvability of the spin chain 



In this section we shall see that the spin chain ([5]) and the spin dynami- 
cal model are related via a slight modification of the usual freezing trick 
mechanism. We shall exploit this connection to compute in closed form some 
eigenstates of 7i with integer energy, for any number of particles and arbitrary 
spin. We shall then compare our results with those obtained by numerical di- 
agonalization of 7i for small values of and n. The numerical computations 
strongly suggest that the exact states derived in this section exhaust all the 
eigenstates of the chain ([5]) with integer energy. 

We shall use in what follows the decomposition 

H = H,^ + 2ah{^) , (28) 

where 

/i(x) = Y^ixi - x,+i)-' (1 - S^,i+i) . (29) 

i 

It shall also be convenient to write the scalar Hamiltonian ([H]) as 

^.c = - E 91 + a^U{^) - a V{^) , (30) 

i 

where the scalar potentials U and V are respectively given by 

f/(x)=r2 + E7 ^ + E7 1 (31) 

■ [Xi — Xi-^-i) ■ [Xi — Xi-i)[Xi — Xi^i) 

^(x) = E 7 ^ ■ (32) 

■ [Xi — Xj+i ) 

From Eq. (l30l) it follows that for large a the low-lying eigenfunctions of ifgc 
concentrate at absolute minima of the potential U in the domain (ITOl) and 
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their energies satisfy [15P6] 



E = a^minf/(x) + 0(a) . 

By comparison with the exact results for the ground state ([9]), one concludes 
that ^ is an absolute minimum of U and U{^) = 2N. This condition has been 
used in the previous section to check the accuracy of the numerical solution 
of the sites equations for large values of A^. 

There are two main limitations which prevent the use of the standard freezing 
trick argument [2511^ for computing the spectrum of the spin chain The 
first one is the requirement that ^ be the unique minimum of the potential 
U in the domain C. Although our numerical calculations suggest that this 
is indeed the case, we have not been able to provide a rigorous proof of this 
fact. The second limitation, which is more fundamental, is the fact that the 
dynamical models H and ifgc are only quasi-exactly solvable. Let us briefly 
recall the basics of the usual freezing trick method in order to understand why 
the full knowledge of the spectra of H and ifgc is essential for its application. 
Indeed, if the potential U has an unique minimum ^ in C, for sufficiently large 
a all the eigenfunctions of Hsc are sharply peaked around this point. Thus, if 
ip{x.) is an eigenfunction of ifgc with energy Esc and |cr) is an eigenstate of the 
chain Ti with eigenvalue S, for a ^ 1 we have 

/i(x)V^(x)|(t) ~ V^(x)/i(^)|(t) = ^{^)n\a) = S^{x)\a) . 

By Eq. f l28|) . the state ?/'(x)|(t) is then an approximate eigenfunction of H with 
eigenvalue 

E Esc + 2a£ . (33) 

In other words, the Hamiltonian ([7]) is approximately diagonal in a basis of the 
form {'i/'j(x)|(Tj)}, where {ipii'x)} is a basis of eigenfunctions of H^c and 
is a basis of eigenstates of H. Equation cannot be used directly to compute 
the corresponding spectrum of Ti, since it is not clear a priori which pairs of 
eigenvalues of H and Hgc yield an approximate eigenvalue {E — Esc) /{2a) of 
Ti. However, using Eq. fl33l) one can easily express the partition function Z of 
the spin chain Ti in terms of the partition functions Z and Zgc of H and Hsc, 
respectively, via the formula 

Z{T) = hm . (34) 

^ ^ ™Z,e(2aT) ^ ' 

The latter formula, which is the key result behind the standard freezing trick 
approach, cannot be used to compute the spectrum of the spin chain H unless 
the whole spectrum of both H and H^c is known. 

In spite of the above limitations, we shall see in this section that it is still 
possible to compute a number of eigenstates of the spin chain ([5]) from some 



13 



of the families of spin eigenfunctions of the model ([7]) constructed in Ref. [M]. 
More precisely, we will show that certain linear combinations of these eigen- 
functions factorize as the product of the ground state of the scalar model ([S]) 
times a spin function, whose limit as a — > oo is an eigenstate of the spin chain 
Hamiltonian H.. 

Let us begin by introducing some preliminary notation. Let S be the space 
of internal degrees of freedom of particles with su(?t,) spin, and denote 
the elements of its canonical basis by |si...SAr), where Si = —M,—M + 
1, . . . , M and M = {n — l)/2. Let A be the total symmetrizer under particle 
permutations, that is 

1 m 

where denotes a permutation operator acting simultaneously on spatial 
coordinates and spins. We shall consider in what follows the subspace S' C S 
of spin states |s) such that J2i is symmetric, where \sij) is defined by 

A(^xiX2\s)'j =^XiXj\sij) ] (35) 

i<j 

see Ref. [3^ for a complete characterization of this subspace. Finally, given a 
state |s) G E we define the spin functions 

$('=)(x;|s)) = A(x^|s)) (A; = 0,1, 2), $(2)(x;|s)) = A(xiX2|s)) , (36) 
$(=^)(x;|s)) = A(a;iX2(xi - X2)\s)) . (37) 

We shall suppress one or both of the arguments of the above spin functions 
when appropriate. 

We shall next present the basic result we have used to construct eigenstates 
of the chain out of the eigenfunctions of the dynamical spin model ([7j) in 
Appendix [Al Recall that the energies of these eigenfunctions are the numbers 
Elm = Eq + 2a(2/ + m), where / and m are non- negative integers and Eq is 
the ground state energy of the scalar model ifgc- 

The key point in the ensuing argument is the fact that for large a the normal- 
ized ground state /xq = /^/||/^|| of Hgc, where is given by and = /i^, 
is sharply peaked around its maximum ^. Hence, if E{x.) is a continuous spin- 
dependent function such that the integral of fi^F over C is finite, the main 
contribution to this integral comes from a small ball centered at ^, up to expo- 
nentially small terms. By the standard argument behind the proof of Laplace's 
method, it follows that 

Jl-i^oi^ = ^(^)il-X/^o = i^(0. (38) 
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Assume now that \I/(x;a) = /io(x; a)$(x; a) is a continuous eigenfunction of 

k 

H with energy Eki such that $(x; a) = ^ a~*$j(x). Denoting by 



j=0 



(/,F)= / /(x)F(x)d^x 



the usual spin- valued inner product of a scalar function /(x) with a spin- valued 
function -F(x), from Eq. (128|) we obtain 

£;m(/^o. = ^fcKAio, ^) = (yuo, = (yuo, ^sc*) + 2a(/io, Z^^) • (39) 

Since ifgc is self-adjoint and /i is a matrix multiplication operator, the RHS of 
the previous equation equals 

(i/sc/io, /io$) + 2a(/i2, = Eo(/io, /io'^') + 2a(/i2, = (/i^, (^^ + 2a/i)$) . 
From Eq. fl39|l it follows that 

= (/i^, (/i - 2/ - m)<l>) = ^ a-*(/i2, (/i - 2/ - m)<l>i) 

Taking the limit a — oo, and using Eq. (138|) and the fact that = h{$), we 
finally obtain 

(7^-2/-m)<l>o(0 = 0. (40) 

Thus 

<l>o(0=ilm<f(0 (41) 
is either zero or an eigenstate of the chain ([5]) with integer energy 21 + m. 

We shall now see that the application of the previous method to the eigen- 
functions flA.ip - flA.SP of the dynamical model ([7]) listed in Appendix \K\ yields 
the following types of eigenstates of the chain ([5]): 

= ^^'\\s)) ^ k\s) , (42a) 
^ii\s)) = . (42b) 

^2(k)) = + {N- , (42c) 

V,{\s)) = ¥^\t\s)) + 2^^^\i-\s)). (42d) 

For the states fl42cl) |s) must belong to the subspace S' and be symmetric in the 
first two spins, whereas for the states of type fl42dp \s) should be antisymmetric 
under 5*12. In all cases, the energy of the states ipk{\s)) is the integer k. We 
shall see in what follows that the states fl42dp of energy 3 only appear for spin 
M > 1/2 {n > 2). 

Note, first of all, that Eq. flA.Sp implies that = A^, since otherwise each 
eigenfunction ^'[q'* (cf. Eq. (lA.ip ) would yield an eigenstate of the chain ([S]) 



15 



with eigenvalue 21 for arbitrary /. (The fact that = can also be es- 
tabhshed directly from Eq. using the identity U{^) = 2N for the scalar 
potential fl?I]) .) Hence we need only consider the eigenfunctions flA.ip - flA.5p 
with / = 0. Moreover, from Eqs. f]A.l|) - f]A.5P and flA.9P it follows that the 



eigenfunctions ^E^om; ^om ^^id \E'o^ (respectively \E'o^ and '^'om) vanish identi- 
cally at X = ^ for odd (respectively even) m. 

For the first type of eigenfunctions f lA.ip . when m is even Eq. f lA.lOp implies 
that 

AT — 

CPomW - a ^flo Wom ^ ^rny/^ ' 

where is defined in (jSS]). By Eq. ([T2D, the RHS of the previous equation 
vanishes at x = ^ unless m = 0. Thus the eigenfunctions flA.ip only yield the 
zero energy eigenstates f l42ap of Ti. (when I = m = 0). Note that the fact that 
any symmetric spin state is an eigenstate of the chain Ti. with zero energy 
follows directly from Eq. fJS]). Conversely, the latter equation implies that the 
ground state energy of Ti. is zero and that the corresponding eigenstates are 
totally symmetric. For this reason we shall concentrate in what follows on the 
nontrivial states fl42bl) - fl42dl) . 



Let us next examine the eigenfunctions of type flA.2p for / = and odd m ^ 1. 
In this case we have 



-r rn~l 

N— 



$i^i(x) . a^^o^^i^ — 7^^"^-^*^^^ ' (43) 



where is given by fl5^ . Using again Eq. flT^ we conclude that lima^oo ^'om 
vanishes at x = ^ unless m = 1. In this case, we have 

lim$«(0 = <l.«(^;|.))^^i(|.)). 



(2) ^ f'2) 

Consider now the eigenfunctions \E'om and \t'o,^ with even m. Using again 
Eq. flA.lOp we have 



2N-X''' 
im-T)(%^)! 



a-^/^o = --^^^7^^ + O(a-) , (44a) 



""^ (iV-l)(m-l)(n^^)! ^ ^ ^ ^ 

where m ^ 2. Thus the 0(1) term of the left-hand sides of Eqs. f j44l) vanish 
at X = ^ for all m ^ 2, on account of Eq. ffT2|) . However, if '^^J^ and '^q^ are 
built from the same spin state \s), the previous equations imply that the 0(1) 
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m 1 / {'?\ 

part of the linear combination a'^jji^ 
observation suggests considering the function 

$(2) + (AT _ i)$(2) _ iVx(2$(^) 



l)^om) vanishes. This 



m-2 



where m ^ 2 is even and ^'^^^ are given in Eq. fl36|) . Note that is only 
defined for states |s) G S' symmetric under S'12, since otherwise would 
not be defined. The limit of $om as a ^ 00 vanishes at x = ^ unless m = 2, 
in which case we obtain the eigenstate with energy 2 



Let us finally turn to the last type of eigenfunctions flA.Sp with / = and odd 
m ^ 3. From Eq. flA.lOp it immediately follows that 

a^/^o-^« = /"^^frl { - - ^''') + Oia-^) , (45) 



(m-2)(^)! V m 



whose limit as a ^ 00 vanishes identically at x = ^ since m ^ 3 in this case. 
However, as for the previous states, one can cancel the leading term in a of 
\l/om with a suitable linear combination of \l/o^ and \&om all built from the 
same spin state (necessarily antisymmetric under S'12 for $'0^ to be defined). 
Indeed, using Eqs. P3|) and fHSl) and taking into account that for odd m 



a^/io^^i°i = ^^-$(°^+0(a-), 
V 2 j- 

one can easily check that 

= a^A„-'(*S. - - S-S) = 0(1) . (46) 

V m — 2 m[m — 2) J 

We shall next see that the 0(1) term of $o^(x) vanishes at x = ^ unless 
m = 3. Indeed, Eq. (IA.9P implies that for k' > 2k the term x^' Pl^^'^'^\t) 
vanishes at x = ^ to all orders in a. Hence the 0(1) terms of $o^(x) and 



a'^x--' ( Ptf''\t) $(3) - 2^ x'Pt\'''\t) $(1)) (47) 
V m — 2 / 

coincide at x = ^. From Eq. (]A.9P it is straightforward to show that the 0(1) 
term of (1471) is a linear combination of x"^~^ and x"^~^, thus establishing our 
claim. 
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By the previous remarks, we need only compute the 0(1) term of (H7|) for 
m = 3, which by Eq. (1A.10|) is given by 



Using Eq. f|T2l) and the identity = we finally obtain the eigenstates fl42dl) 
of energy 3: 



Jim <|.g)(0 = ^^'\$;\s)) + 2$«(^;|.)) ^ v,3(|^)) • 



The previous discussion guarantees that the nontrivial states ipk{\s)) given by 
Eqs. (I42bl) - fl42dl) are eigenstates of the chain H with energy k provided that 



they do not vanish. For instance, the states fl42bp are easily seen to vanish 



when the spin state \s) is symmetric, since in this case 

This is also the case for the states fl42cp . Indeed, if \s) is symmetric it clearly 
belongs to S' and we have 

$(^)(x;|s)) + (iV - l)¥'\^;\s)) = ^ |s) + ^ \s) = Nx'\s) . 

i<j 

Since ip2{\s)) coincides with the LHS of this expression evaluated at x = ^, it 
vanishes on account of Eq. (fT2l) . Less trivially, let \s) be a linear combination 
of basic states |s) = |si . . . sjy) such that |s| ^ 2, where |s| denotes the number 
of distinct components of s = (si, . . . , sat). If, in addition, |s) is antisymmetric 
under Su, we shall prove below that 



2 

$(3)+2(^<l>«-x<l>(2)) =0, (4J 



which implies that (p-i{\s)) = by Eq. ( fT2|) and the identity = A^. In 
summary, the state \s) in the definition of the states ( 142bl) -( l42dl) can be taken 
without loss of generality as follows: 

^i(|s)): Is)GS-A(E); (49) 

Is)gS'-A(S), Su\s) = \s); (50) 

^3(k)): |s)=^c,|s), Si2\s) = -\s). (51) 

|s|^3 



Note, in particular, that the states (142dp do not appear for spin 1/2 in view 
of the first condition (15T!) . 

Our next task is to study the number of linearly independent states of the 
form (142bp - (142dp with \s) satisfying the above conditions. The main difficulty 
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in this respect is the fact that the eigenstates ipk{\s)) constructed from a set of 
hnearly independent states |s) satisfying conditions (149|) - (!5T]) need not be in- 
dependent. In order to address this problem, it is convenient to introduce some 
additional notation. Given a basic state |s) with |s| = p, we define its spin 
content as the set of pairs (s*, Ui), where < ■ ■ ■ < are the distinct compo- 
nents of s and z/j is the number of times that appears in s. For instance, the 
spin content of the basic state |s) = |0, —2, 1, —2, 1) is {(—2, 2), (0, 1), (1, 2)}. 
We shall say that an arbitrary spin state \s) has well-defined spin content if 
it is a linear combination of basic states having the same spin content. It is 
obvious that a set of states whose spin contents are all different is linearly 
independent. Note also that if \s) has a well-defined spin content, then ipk{\s)) 
has the same spin content. Therefore, it suffices to determine the number of 
linearly independent states ipk{\s)) built from states \s) with well-defined spin 
content. Clearly, for a given spin content {(s*, z^j)} this number is independent 
of the particular values of the spin coordinates s*. We shall therefore denote 
by (ifc(z/i, . . . , i/p) = dki^v) the dimension of the linear space of states ipk{\s)) 
with a given spin content {(s*,z/j)}. We shall next prove the following upper 
bounds on these dimensions: 

diH^p-l, daM^p-l, 4M^r2 j- ^^^^ 



Consider in the first place the states of the form (fi{\s)) with a given spin 
content. If |s) and |s') are two basic states with the same spin content and 
differing by a permutation of the last — 1 spin coordinates, then A(a;i|s)) = 
A(xi|s')), which implies that v^ids)) = (pi{\s')) by Eqs. fl36|) and fl42b|) . Hence, 
the space of states ipi{\s)) with spin content {{s^ , ui) , . . . , (s^, z/p)} is spanned 
by the states 

M\s\..)), z = l,...,p, (53) 
where the ellipsis denotes any ordering of the remaining A^— 1 spin components 
corresponding to the above spin content. On the other hand, the p states (1351) 
satisfy the linear relation 

j2i^iM\s'...)) = 0, (54) 

i=l 

which implies the first inequality in (!52|) . Indeed, if |s) is a basic state with 
spin content {{s\ i^i)}, we have 

N p 

NxA\s) = J2a{x,\s)) = j2 ^^^{xl\s' •••))• 

i=l i=l 

Setting X = ^ and using (fT2l) we immediately obtain (15^ . 

Let us turn now to the states of the form ip2{\s)), where \s) satisfies ( !50!) 
and has a well-defined spin content {(s^, z/i), . . . , (s^, Pp)}. The results of our 
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previous paper [M] (see Proposition 4) imply that the space of states of this 
form is spanned by 

^2(|Xi)), ^ = 1,---,P, (55) 

where 



\Xi) 



(56) 

z/,(|sV...) + |sV...)), u, = l. 



However, from jMi Prop. 3] it follows that the p states (ISSj) satisfy the linear 
relation 

^2M\X^))=0, (57) 

i=l 

which yields the second inequality in (!52l) . 

Consider finally the states of the form ipsds)), where \s) is antisymmetric 
under 5'i2 and has a certain spin content {(s^, ui), . . . , (s^, t'p)}. Clearly, the 
space of such states is spanned by 



^3(|sV...)-|s^s\..)j, l^t<j^p. (58) 
We shall prove below that the latter states satisfy the linear relations 

f]z/,<^3(kV...)-|s^s\..)) =0, z = l,...,p. (59) 

Multiplying the LHS of the i-th relation by Ui and summing over i we obtain 

p 



which vanishes by antisymmetry. Thus we can safely drop one of the iden- 
tities dsn]), say the last one, and the first p — 1 are independent since they 
can be solved for ips(^\s^s^ ■ ■ ■) ~ ■ ■ ■ )) i = 1, . . . ,p — 1. Hence there 
are at most — {p — 1) = (^2^ independent states of the form fl58|) . which 
establishes the last inequality in fl52|) . 



We still have to prove the identities fl59|) . To this end, denote by $3 = $3(x; |s)) 
the LHS of Eq. (HH}. Since $3(^;|s)) = (^3(|s)), Eqs. follow directly from 



the identities 

p 



J2 ^z\^¥s^ ■■■)- \s's' ...))= , z = 1, . . . ,p, (60) 
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In order to establish (160|) . note first that 

N 

r2$(i)(|s^s^..)) - Nx^^^^{\s's^ ...)) = - ^A(xiXi(xi 

1=1 

= -t.,$(^)(|.V...)), 
fc=i 

so that 

^ ^2 ^3{\s's^ ...)- \s^s' ...))= Nj^i'j ^^'^ {\s's' ...)) 

- X: ^ z/,z/,l>(3)(|s^s^..)). (61) 

j,A:=l j,k=l 

Since Z]j=i = the first two terms of the RHS of the previous formula 
cancel, whereas the last one vanishes by antisymmetry. This concludes the 
proof of Eq. fl60l) . Note, finally, that the latter equation for p = 2 easily yields 
the identity fHg]l . 

3. 1 Example 

We shall now discuss in detail the case = 5 and spin M = 1. The sites 
of this chain are given by ^5 = —,^1 = 1.469, ,^4 = —^2 = 0.584, ^3 = 0. The 
Hamiltonian ([5]) is represented by a 243 x 243 real symmetric matrix, whose 
numerical diagonalization is straightforward. The spectrum consists of 21 dif- 
ferent levels, the only integer energies being 0,1,2,3. For conciseness' sake, 
we only present the first six levels and the last one with their corresponding 
degeneracies, cf. Tabled Notice, in particular, the appearance of a noninteger 
energy between the third and the fifth levels. The degeneracy of the ground 
level is given by the number of independent symmetric states, that is, the 
number of combinations with repetitions of 5 elements from 3. Turning to the 
other integer energies, we shall next present a basis of eigenstates of each non- 
trivial type fl42bp - fl42dp . By the previous discussion, we must first determine 
all the possible spin contents {(s*,z/j)} compatible with conditions P9l) - fl5T|) . 
and then construct a basis of states ipk{\s)) with each of these spin contents. 
In practice, the former task is performed in two steps, namely, one first finds 
all possible degeneracy vectors u, and then one determines the spin contents 
associated with each of these vectors. 

The states (142b|) of unit energy are generated by basic states |s) with p ^ 2 
different spin components. In Table [2] we present the list of such basic states, 
where we have taken into account the relation (15^ to drop one state for each 
different spin content; for instance, the state |0 ) does not appear in the 
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Table 1 

Energies and degeneracies 
iV = 5 and M = 1. 



of the first six and the last levels of the chain ([5]) with 



Energy Degeneracy 






21 


1 


24 


2 


24 


2.035 


15 


3 


6 


3.953 


15 




15.033 


3 



table, since v^idO )) = — 4(/?i(| 0)). In this way we obtain 24 basic 

states, whose corresponding states {pi{\s)) turn out to be linearly independent. 

Table 2 

List of basic states |s) generating the states ()42bp of unit energy for each degeneracy 
vector u. 



V 








(4,1) 


0), 


— +), 


I0000+) 


(3,2) 


00), 


— ++), 


I000++) 


(2,3) 


000), 


— +++), 


I00+++) 


(1,4) 


-0000), 


-++++), 


I0++++) 


(3,1,1) 


0+), 


+) 


(2,2,1) 


— 00+), 


00 +) 


(2,1,2) 


— 0++), 


++) 


(1,3,1) 


-000+), 


000-+) 


(1,2,2) 


-00++), 


00-++) 


(1,1,3) 


-0+++), 


0-+++) 



Similarly, the space of states V52(|s)) of energy 2 is spanned by the states ( l55i) 
with \xii given by ( l56i) . for each spin content {(s^, z/i), . . . , (s^, z/p)} such that 
j9 ^ 2. As in the previous case, the relation (1571) implies that for each spin con- 
tent one can drop one of the states In Table[3]we present a possible choice 
of the 24 states that can be constructed by this procedure. Just as before, 
the corresponding states v^2(|Xi)) are easily seen to be linearly independent. 
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Table 3 

States I x) of the form ([56]) generating the states (j55|) of energy 2 for each degeneracy 
vector V. 



u 




Ix) 


(4,1) 


1 o),| 


+), I0000+) 


(3,2) 


3| 00) - 

3|000++) - 


-|00 ),3| ++)-!++ ), 

I++000) 


(2,3) 


3|000 ) 


1 000) , 3|+++ ) - 1 +++), 


3I+++00) - 


- I00+++) 


(1,4) 


I0000-), I++++0) 


(3,1,1) 


6| 0+) - 

3(1 0—+) 


_|0+ )-| + ), 

+ 1-0 — +)) + |o+ ) + | + o ) 


(2,2,1) 


1 00+) - 

|00 +) - 


|00 +) - |0+ 0) - 1 + 00 ), 

1 — 00+) - I-+-00) - 1+ 00) 


(2,1,2) 


1 0++) - 

I++ 0) - 


I++ 0) - |0++ ) - 1+0 +) , 

1 ++0) - I-0-++) - |0 ++) 


(1,3,1) 


3(1-000+) + I0-00+)) + I-+000) + I+-000) , 
6|000-+) - I-+000) - I+-000) 


(1,2,2) 


I00-++) - 
I++-00) - 


I++-00) - I-++00) - I+-00+) , 
I00-++) - I-00++) - I0-0++) 


(1,1,3) 


I-0+++) + 
6I+++-0) - 


I0-+++) +3(|-+++0) + I+-0++)) , 
- I-0+++) - I0-+++) 



Finally, for spin 1 the set of states ipz{\s)) of energy 3 is spanned by states of 
the form (!58|) associated with spin contents with three different components, 
cf. Eq. ( [51]) . By the last inequality (!52l) . for each such spin content there is 
at most one independent state of the form (1581) . Thus there are at most 6 
independent states of the form (]42dl) . generated (for instance) by the states 
\s) listed in Table HI As in the previous cases, it can be checked that the 6 
states ip^(\s)) constructed from the states in Table [Hare actually independent. 

Several important remarks can be made in connection with the previous ex- 
ample. In the first place, it is apparent that the inequalities fl52l) are in this 
case equalities for every spin content. Secondly, the number of independent 
states with integer energy k of the form yj^ds)) coincides with the degeneracy 
of the corresponding level, cf. Table [1] Hence, in this example all the eigen- 
states with integer energy are of the form (H2|) . and can therefore be computed 
explicitly. In fact, we have performed a similar study for the case of = 6 
and spin M = 3/2, arriving at exactly the same conclusions. In view of these 
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Table 4 

List of states |s) generating the states ()42dp of energy 3 and their corresponding 
degeneracy vectors u. 



V 



(3,1,1) 


|0+— ) 


-1 + 0—) 


(2,2,1) 


1+ 00) 


- I-+-00) 


(2,1,2) 


|0 ++) 


- I-0-++) 


(1,3,1) 


I-+000) 


- I+-000) 


(1,2,2) 


I-00++) 


- I0-0++) 


(1,1,3) 


I-0+++) 


- I0-+++) 



examples, it is natural to formulate the following conjectures: 

1. The inequalities ( l52l) are always equalities. 

2. The only integer energies of the chain ([5]) are 0,1,2 and (for M ^ 1) 3, 
corresponding to the first three and the fifth levels. 

3. The only eigenstates with integer energy are those of the form (l42l) . 

According to the first conjecture, for = 1, 2, 3 the number of independent 
states 9?fc(|s)) with a well-defined spin content with p elements depends only 
on p, and is given by 



4(p) 



>-l. A; = 1,2 



2 ) ' 



k 



By the third conjecture, the degeneracy of the integer energy k is thus 

A; = 1,2, 3. (62) 



mm{n,7V) (n\(N-l' 

dk= Yl dk{p)\ 
p=i 



.pj \P — 1 



Indeed, there are (^^^ different choices of p distinct spin values s*, and (^Ji 
ways of selecting p numbers z/j ^ 1 such that ui + ■ ■ ■ + Up = N . In order to 
evaluate the sum in (l62l) . note first that 



min(?i,A'') / iX/X/at i^ 
' I p — 1\ I n\ I N — 1 

p=l 



p=i 



j )\p)\p-i) \ J J \pj\p-j-i 

On the other hand, expanding both sides of the identity 



(63) 
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in powers of z, we easily obtain 



minfrtjA'") / \ /at ■ ^\ 



E 

p=i 



.pj\p-j-i 




The previous identity and Eq. fl63|) immediately yield the following explicit 
formula for the degeneracy of the positive integer energies of the chain ([5]): 

(The degeneracy of the ground level £^ = is the dimension of the space of 
symmetric states, namely (^^~^^~^^ •) 

We have numerically diagonalized the chain Hamiltonian ([5]) for several values 
of n and N such that n'^ ^ 3^ = 6561. In all cases, we have checked the 
validity of the second conjecture, and that the formula flMl) for the degeneracies 
of the integer energies (which is a direct consequence of the first and third 
conjectures) is satisfied. These results lend strong numerical support to the 
above conjectures, whose rigorous proof deserves further investigation. 



4 Statistical analysis of the spectrum 



It has been recently shown that for large N the level density of spin chains 
of HS type can be approximated with remarkable accuracy by a Gaussian 
distribution [29|31|32j . A natural question is whether this is also the case for 
the spin chain under consideration. The fact that the partition function of 
the chain has not been computed in closed form — unlike those of the HS 
chains in the previous references — makes it difficult to address this problem 
directly. Although one can in principle diagonalize numerically the matrix of 
the chain Hamiltonian Ti, in practice this is only feasible for relatively small 
values of A^. In this section we shall study the level density of the chain ([5]) 
using the method of moments outlined in Appendix [Bl which provides an 
accurate estimation of the level density for larger values of A^. 

In order to improve the numerical stability of the method, it is convenient to 
take the parameters b and c in Eq. flB.ll|) as 

b = 0, c = 1 , (65) 

so that the spectrum of A lies in the interval [—2, 2]. Therefore the matrix A 
and the Hamiltonian Ti are related by 

AH 

A = - 2, (66) 

^max 
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where S^ax is the largest eigenvalue of H. Note that this eigenvalue can be 
computed numerically without difficulty for relatively large values of A^, since 
the matrix of Ti is very sparse. One must then check that the coefficients bk and 
Cfc in the continued fraction expansion of the resolvent of A (cf. Appendix [B|) 
approximately stabilize to the values b and c in Eq. fl65l) for ko < k < ki, with 
ki ^ fco- We have verified that this is indeed the case for sufficiently large, 
where ko is typically of the order of 10. We can thus use Eq. flB.12|) to obtain 
a continuous approximation g{x), where x = {4S/Smax) — 2, to the density of 
eigenvalues of the matrix A, as explained in Appendix [Bl 

As a test of the accuracy of the method, we have computed the density g for 
spin 1/2 and A^ = 12, since in this case the matrix A can still be diagonalized 
numerically. In Fig. O we have compared the approximate density g computed 
by applying the moments method with ko = 20 and 20 random vectors with the 
histogram of the spectrum of A obtained by subdividing the interval [—2, 2] in 
50 subintervals. As can be seen from the latter figure, the continuous density 
g essentially reproduces the shape of the histogram. 
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Fig. 6. Continuous approximation g{x) to the density of eigenvalues of the ma- 
trix (|66|) for spin 1/2 and = 12 {E-max = 201.21) compared with the histogram 
of the spectrum of this matrix. 



The real interest of the method of moments is the possibility of approximating 
the density of eigenvalues of a large matrix whose numerical diagonalization 
is not feasible. We have been able to compute the approximate density g{x) 
of eigenvalues of the matrix A for A^ up to 21 for spin 1/2, and up to 13 
for spin 1. For instance, in Fig. [7] we present the plot of g for spin 1/2 and 
A^ = 21, computed with ko = 20 and 20 random vectors. It is apparent from 
this plot, and for similar plots for spins 1/2 and 1, that the level density is 
not Gaussian. As a matter of fact, the function g is well approximated by the 
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Wigner-like distribution 



/(x)=Ar-Ve-^, y = x + 2, 



(67) 



where the normahzation constant M is given by 

Af = -„,3.-,r(-),F.(— .-;-r)+_^^r(— 

and iFi is the confluent hypergeometric function of the first kind. Thus the 
behavior of the chain Ti is rather different in this respect from that of the spin 
chains of HS type. 

f{x),g{x) 




\ 



J 



Fig. 7. Approximate density g(x) of eigenvalues of tfie matrix (j66p (grey line) for 
spin 1/2 and = 21 (<?max = 1126.53), compared with the distribution ([67|) with 
optimal parameters a = 0.85, (3 = 1.27, 7 = 0.91 (dashed line). 



When N is greater than 12, the accuracy of the level density obtained by the 
method of moments cannot be gauged directly by computing the spectrum of 
the matrix A. An indirect way of estimating this accuracy consists in compar- 
ing the first few moments of the density g with those of the spectrum of A, 
computed by taking the traces of appropriate powers of Ti. As an example, we 
shall next derive simple expressions for the mean and variance of the spectrum 
of Ti in terms of finite sums involving only the coordinates of the chain sites, 
which can be easily evaluated for very large values of A^. We shall see in this 
way that the agreement between these values and those obtained from g is 
indeed very good and, roughly speaking, improves as A^ grows. 

In the first place, the mean energy £ oiTi can be easily computed noting that 
tr5i,i+i = so that (cf. Eq. ([22])) 

E = n-''tin = Hn"" - n""'^) = (l - ^) E • (68) 
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Similarly, since 



and 



n 



N-2+25.. 



'J 



a straightforward calculation yields 

tT{n') = {n- l)n^-^ f (n - 1) ( E /^O' + + 1) E 



so that the variance of the energy is given by 



a^^n-^tr(7^^)-^^=(l-l)E/^^ (69) 



In Fig. [8] we present a logarithmic plot of the relative errors 

between the exact values fl68l) - fl69l) and their approximations 

£, = ^ j\ + > = % /_V + 2)'^7(x) dx - E] (71) 

for = 12, ... , 21 and spin 1/2. From the latter plot it is apparent that for 
N ^ lb both errors are less than .1%, which suggests that for large N the 
continuous function g[x) computed by the method of moments is indeed an 
excellent approximation to the level density of the matrix A. 

From Eqs. fl68|) and fl69|) and the discussion of the distribution of the chain 
sites in Section [2] it is straightforward to deduce the asymptotic behavior of 
the mean and the variance of the energy for large N . Indeed, using Eq. f l24l) 
we easily obtain 

hu ^ fe+i - ik)-' ^ [x{k + 1) - x{k)]-' ^i^—) ^ — e-?i , (72) 
so that 

_ ,2 ^ ,^ ^ ^ N _il 

2 



k 27r V 2vr Y 



e 2-^ dx 



(27r)3/2 7-00 27rv^ ' 
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15 16 



18 19 20 21 



N 



Fig. 8. Logarithmic plot of the relative errors Ai (black boxes) and A2 (crosses) 
in Eq. (|70p as a function of the number of sites N for spin 1/2. The approximate 
level density g in Eq. (|71|) has been computed taking = 20 and averaging over 
20 random vectors. 



Hence for large N the mean energy S behaves as 

Ar3 



04) 



(73) 



Similarly, 



k k 



N 



— Ik 

e 2 



ATS 



e 2-^ dx 



ATS 



^27r (27r)5/2 47r2v^' 

and therefore the large N limit of the variance of the energy is given by 



ATS 



47r2v^ 



1-^1. 



(74) 



Thus, the leading behavior of the energy mean and variance of the chain is 
analogous to that of the trigonometric HS spin chains studied in Refs. [2H1I3I]- 
In particular, both quantities grow with much faster than their counterparts 
for the Heisenberg chain (both of which diverge as A^). 



5 Conclusions 



In this paper we have studied a novel spin chain with position-dependent near- 
est neighbors interactions, which is intermediate between the Heisenberg chain 
(position- independent, nearest neighbors interactions) and the spin chains of 
Haldane-Shastry type (position-dependent, long-range interactions); see Ta- 
ble O for a brief comparison. 
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We have developed a new method related to Polychronakos's freezing trick 
which has made it possible to compute in closed form a certain number of 
energy levels with their corresponding eigenstates, for any number of particles 
and arbitrary spin. While the eigenvalues of the original Haldane-Shastry 
and Polychronakos-Frahm chains (the closest HS analogs of the chain under 
consideration) are all integers, the only integer eigenvalues of our chain are 
precisely those that have been computed exactly by our method. This fact 
strongly suggests that the chain ([5]) should be regarded as the first instance of 
a quasi-exactly solvable spin chain, thus extending the usual notion of quasi- 
exact solvability to finite-dimensional Hamiltonians. 

Table 5 

A comparison of some properties of the spin chain TC versus the Heisenberg chain 
TCuc and the Polychronakos-Frahm chain TCpp (the properties followed by an aster- 
isk are based on unpublished work by the authors). 





n 




TCpp 


Chain sites distribution 


Gaussian 


Equispaced 


Circular law (I26p 


Solvability 


Quasi-exact 


Exact 


Exact 


Integer energies 


First few 


No 


All 


Level density 


Wigner-like law ()67p 


Wigner-like* 


Gaussian* 


Mean energy growth 




N 


Ar2* 


Energy variance growth 


iV5 


N 


Ar3* 



The method developed in this paper is quite general, and only relies on the 
explicit knowledge of a number of eigenfunctions of the corresponding spin 
dynamical model. In this respect, it goes one step beyond the usual freezing 
trick, which requires the computation in closed form of the whole spectrum of 
the related dynamical model. In particular, our method is well-suited to spin 
chains whose associated dynamical model is only quasi-exactly solvable, like 
for instance the models with elliptic interactions constructed in Refs. [1711^ . 
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A Exact eigenfunctions of the dynamical spin model ([7]) 



In this appendix we shall list the eigenfunctions of the dynamical spin model ([7]) 
used to construct the eigenstates of the chain presented in Section [31 Let 



a = N{a + 



3 
2 



p = (3{m) = l-m- N[a + 



2r^ 



- 1 



where x = -^^iXi is the center of mass coordinate. In Ref. [M] it was shown 
that the model ([7j) possesses the following families of spin eigenfunctions with 
energy Eim = Eq + 2a(2l + m) , where / ^ and m is as indicated in each case: 



x™Lr''(ar2)PfSf (t) , m ^ , 



m ^ 1 



iTf(2) _ ,, ™m-:^ r 



-m-2 r -13, 



ar 



2(a + i: 



2[^J + 



m > 2 



(A.l) 
(A.2) 



(A.3) 



+ 



2(a + 1) 



(2L^J+l)(A^-l) 



m ^ 2 , (A.4) 



iT/(3) _ ™m~3 r -/3, 



ar 



L 2 J ^ 3 



2(2a + 3) 



(0) 



mim 



m — 3 
2 



,(0) 



m ^ 3 - 



(A.5) 



In the previous formulas [-J denotes the integer part and e(m) = |(1 — (—1)'") 
is the parity of m. The spin functions ^^'^^ $(2), $(3) (cf. Eqs. (EE]) and ([371)) 
are built from a state |s) symmetric under 5*12 and belonging to S' for "^i^, and 
antisymmetric under Su for The generalized Laguerre polynomials L'[^ 
and the Jacobi polynomials P^^'^"^ appearing in Eqs. (IA.ip - (IA.5P are defined 
as 



-Pi 



3 j! 



(A.6) 



P 



(7,/3), 



u E ^ (-^).(7 + /3 + A: + 1),(7 + J + l)fc-, (1 - ^)^' , (A.7) 

• 7=0 -I' 
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where 

{x)j = x{x + 1) ■ ■ ■ (x + j — 1) 

is the Pochhammer symbol. The eigenfunctions of type flA.51) are independent 
of the remaining ones flA.ll) - flA.4p only for spin greater than 1/2 (n > 2), 
cf . [21] . The model ([7]) possesses two additional families of eigenfunctions (also 
derived in the previous reference), which have not been listed above as they 
do not yield any further eigenstates of the chain ([5]). 

Let us now discuss the behavior of the terms L~[^{ar'^) and x"^^ pjf~^^'^\t) in 
Eqs. f lA.ip - flA.Sp as a ^ cxD. Consider first the polynomial LJ^{ar'^). Since 

it follows that 

a-'Li\ar^) = ^Jl^pL + 0{a-^) , (A.8) 

where the term 0{a~^) in the previous equation is actually a polynomial in 
a~^. (Throughout this paper, the symbol 0(a"'^) denotes any function /(a) 
such that f{a) has a finite (possibly zero) limit as a — oo.) For the other 
type of term, note that 



is clearly a polynomial in x. Taking into account that 

{a + i+j + l)fc„, = {Na^'^ (l + O(a-i)) , 

we obtain 



(A.9) 



a X 



k^2kp(aM,0)^^-^ = + O(a-i) , (A.IO) 



k 



k\ 

where 0{a^^) is polynomial in and x 



B The method of moments 



The method of moments [39P0] is a powerful tool for computing the density of 
eigenvalues of a large Hermitian matrix whose spectrum is not known explic- 
itly. The method is based on the relation between a probability distribution 
g{x), i.e., a nonnegative function whose integral over its support [xi, 0:2] is one. 
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and its resolvent 



R(z) = r dx, Imz^O. 

J x\ Z — X 



Using the well-known identity 



lim — — — - = 7i6{x) , 



it is straightforward to show that 

g{x) = ^— lim ± ie) , xE{xi,X2). 
vr e->o+ 

The previous formula makes it possible to compute the probability density 
g{x) if the resolvent R{z) is known. One of the key ingredients of the method 
is the fact that if g is positive in a set of nonzero measure, the resolvent can 
be expanded as a continued fraction 

R{z) = , (B.2) 

z-bo 7 

z ~ bi — ■ ■ ■ 

where bk, Ck are the coefficients in the three-term recursion relation 

Pk+i{x) = ix- bk)Pk{x) - CkPk-i{x) , (P_i = , Po = 1 , Co = 1) (B.3) 

satisfied by the orthogonal polynomial system {Pk(yX)}'^^Q associated with the 
density g{x). The polynomials P^ are the monic polynomials determined by 
the orthogonality condition 

{PkPi) = r Pk{x)Pi{x)g{x) da; = , k^l. 



It can be shown [l9| that these polynomials satisfy a three-term recursion 
relation of the form flB.31) with coefficients b^ and given by 

, {xPi) (Pi) 



In this paper we are interested in computing the density of eigenvalues (nor- 
malized to unity) of a Hermitian d x d matrix A, i.e., 

1 

i=l 
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where Ai ^ . . . ^ are the eigenvalues of A. The corresponding resolvent 

i=l 

is thus a rational function of z. In this case the resolvent can also be expanded 
as a continued fraction of the form (]B.2I) terminating at level d, i.e., q = 
and Eqs. (iRSj) and flRil) hold for A; ^ - 1. 



If the density (7 is not known, the coefficients bk and Ck cannot be computed 
from Eqs. (]B.4p . and hence the expansion (\B.2\i cannot be used directly to 



evaluate the resolvent. On the other hand, the coefficients bk and can be 
determined once the moments (x^) of the distribution g are known. Another 
central idea of the method consists in replacing {x'^) by the average of the 
expectation values (v, ^'^v) over a suitable set of normalized random vec- 
tors V e W^. More precisely, let Ui, i = 1,. . . ,d, be d independent random 
variables uniformly distributed in the interval [—1, 1], and let v be the vector 
with components Vi = Ui/\\u\\. Clearly, the components of v satisfy 

1 



ViVj = - 5ij , (B.5) 

where the bar stand for the average over the random numbers Ui. We thus 
have 

i,j=l i=l i=l 

= i E / ^^^(^ - ^i) = {x'') . (B.6) 
d -^^ J 

The previous equality implies that 



(P(x)) = (v,P(^)v), 



where P is an arbitrary polynomial. From Eq. flB.4p we immediately obtain 
the following formula for the coefficients bk and c^: 



6. = -^^^^^. c. = J^. (B.7) 

llWfclP l|Wfe-lP 

where = Pk{A)\. By Eq. (IB.Sp . the vectors satisfy the recursion relation 

Wfe+i = (A- bk)wk - CkWk-i (w_i = , Wo = V , Co = 1). 

The latter formula, together with Eq. (IB. 71) . can be used to recursively com- 
pute the coefficients bk and Ck for all /c = 1, . . . , ci — 1. 



34 



The number d = is typically very large, so that it is not feasible to com- 
pute all the coefficients and by the procedure just outlined. In practice, 
one only computes the ffist few of these coefficients, say up to k = ko, and 
tries to estimate the part T^oiz) of the continued fraction (lB.2p involving the 
remaining coefficients. In many cases of interest |50j, the coefficients bk and 
Cfc approximately stabilize for a wide range of values of k, namely bk — b and 
Cfc — c for ko < k < ki < d, where ki ^ ko. In this case the remainder Tk^iz) 
can be approximated as follows: 

TtM s T(z) . 

^ 

z — b — ■ ■ ■ 

But the truncation factor T{z) satisfies 

z — b — T[z) 

so that 



T{z) = -[z-b±^{z~bY-Ac). (B.8) 
We thus obtain the following approximate formula for the resolvent: 

R{z) ~ ^- . (B.9) 

z-bQ 

z — bi — 



z - bk„ - T{z) 



The previous equation can be written as 

A{z)±\B{z)f;{z) 
R[z) ^ ^^=, (B.IO) 

C{z)±iD{z)^uj{z) 

with A,B,C,D polynomials in z with real coefficients and 

LU (z) = {z - xi){x2 - z) , xx^2 = b±2\fc\ (B.ll) 

note that x\^2 ^ IR, on account of Eq. ( IB. 71) . The RHS of ( IB.lOp has a branch 
cut on the interval [xi,X2], which is the support of its corresponding den- 
sity function. From Eqs. (IB.ip and (IB.lOp we finally obtain the following ap- 



proximate formula for the normalized density of eigenvalues of the Hermitian 
matrix A: 



\A{x)D{x)-B{x)C{x)\ fAx) 
9{x) ^ \2 ^ -nf \2 f \ ' xe[xuX2]. B.12 
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Since the RHS of the previous equation cannot vanish identically in any subin- 
terval of [a;i,a;2] with nonzero length, it follows that the formula just obtained 
is only appropriate when the spectrum of A has no gaps. 
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